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Abstract 

We explore the correlations between primordial non-Gaussianity and isocurvature per¬ 
turbation. We sketch the generic relation between the bispectrum of the curvature per¬ 
turbation and the cross-correlation power spectrum in the presence of explicit couplings 
between the inflaton and another light field which gives rise to isocurvature perturbation. 
Using a concrete model of a Peccei-Quinn type field with generic gravitational couplings, 
we illustrate explicitly how the primordial bispectrum correlates with the cross-correlation 
power spectrum. Assuming the resulting /nl ~ we find that the form of the corre¬ 

lation depends mostly upon the inflation model but only weakly on the axion parameters, 
even though /nl itself does depend heavily on the axion parameters. 


1 Introduction 


Inflation BE] is widely considered as the best candidate to explain the origin of the primordial 
fluctuations in our universe, and has been well tested by cosmic microwave background (CMB) 
observations, most recently by the Planck mission [3]. The simplest scenario involves a single 
scalar held slowly rolling down a hat potential [2]. However, when embedded in ultraviolet 
complete theories such as string theory |1], it is not likely that such a simple description of 
single held inhation remains legitimate throughout the whole inhationary epoch, as in general 
additional degrees of freedom become dynamically relevant during inhation. 

It is illustrative to consider the held trajectory in the held space to see how these additional 
degrees of freedom ahect the inhationary predictions. While in single held inhation we have 
only a single direction along which the inhaton moves straightly, in multi-held space we have 
more than one direction and thus the trajectory is in general curved. As the component of held 
huctuations along the trajectory is associated with the curvature perturbation, those orthogonal 
to the trajectory are responsible for the isocurvature perturbations [5j. What is important here 
is that since the held trajectory is in general curved, the huctuation along the trajectory at a 
moment receives contributions from those orthogonal to the trajectory before. In other words, 
the curvature and isocurvature perturbations are correlated nm- 

Thus, on general ground the existence of isocurvature perturbations can be interpreted as a 
proof of additional degrees of freedom in the early universe. Even if the isocurvature perturba¬ 
tions themselves are too small to be detected, their correlation to the curvature perturbation 
can lead to distinctive observational signatures. It is thus useful and important to study how ob¬ 
servables are correlated with and sourced by isocurvature perturbations on completely general 
ground. In this article, we concentrate on the correlations between primordial non-Gaussianity 
and isocurvature perturbations due to interactions at horizon-crossing, whereas previous stud¬ 
ies focus on correlations due to non-linear evolution on super-horizon scales [7j. The primordial 
non-Gaussianity we are considering might not necessarily be that associated with the GMB 
scales, but also on smaller scales which are relevant for large scale structure for example. By 
studying the correlation structure, we can put further constraints on the model parameters 
compared to that coming from studying the observables individually. 

This article is organised as follows: in Section [2, we discuss the general structure of three- 
point interaction terms between the inflaton and an isocurvature held, and how the resulting 
bispectrum Bti correlates with the cross-correlation power spectrum Vc- We then illustrate 
this general feature using a concrete example involving a Peccei-Quinn (PQ) type held, which 
interacts with the inhaton via a gravitationally induced coupling in Section [3l We subsequently 
conclude in Section 01 Detailed derivations of the results are given in the Appendix. 

2 General structure of correlations 

For simplicity, let us consider a generic two-held system, one being the inhaton 0 and the 
other being a spectator held 6 that survives after thermalisation and becomes responsible for 
the isocurvature perturbation X at late times. We also assume the curvature perturbation 71 
is purely sourced by the inhaton huctuations. In such a system, the three-point correlation 
function of the curvature perturbation 71 may receive additional contributions from 6 [8] via 
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the quadratic interaction ~ which is responsible for the cross-correlation power spectrum 
Vc- The resulting bispectrum of the curvature perturbation Bji is therefore correlated with Vc, 
giving additional constraints to the spectator held 6. 

To see this explicitly, we decompose the cubic interaction Hamiltonian as 

^ j:^hrossl) jj -{ ctoss 2 ) 

where ~ ^ gi9(l)^, ^hross 2 ) ^ g29‘^(f) and ~ g^9^ with giS being 

the couplings. Note that these are schematic forms only, which could also include derivative 
interactions such as and 9^ . In standard slow-roll inhation usually contributes 

to a slow-roll suppressed, negligible non-Gaussianity of 0{e) [9]. We therefore neglect 
here. Through the transfer vertex H 2 , each cubic interaction Hamiltonian term contributes to 
the primordial bispectrum Bti as 

^ Bn r-g2 V^Vl , ( 2 ) 

^ Bn-- gzVl . 

Diagrammatic representation of each term is shown in Figure [B However, the detail of each 
contribution, snch as the momentum-dependent shape and amplitude, depends on the coupling 
and derivative structure of the corresponding cubic interaction Hamiltonian. In the following 
section using a concrete example we demonstrate that indeed the generic relation ([2]) is valid. 
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Figure 1: Feynman diagrams for the contributions of given by ([T]) to the inflaton three-point 
function. is the intrinsic three-point inflaton self-interaction, which is of 0{e) and is 

usually slow-roll suppressed [9]. ^^cross 2 ) contribute via interactions with 

the isocurvature held. 

Before moving on to discuss a concrete example, let us comment on the relative size of 
the various bispectrum contributions above. Naively, one may think this somewhat suggests 
the above bispectrum contributions are small compared to that from the intrinsic inhaton self¬ 
interaction since 1 -^ Bn ~ Vn, if Vn > Vc- Besides, one might think there exists 

a hierarchy between the various contributions, with However, this 

is not necessarily the case. In general, the resulting contributions are momentum and shape 
dependent. Depending on the forms of and the various contributions 

might all have different momentum dependences and peak at different shapes. Furthermore, 
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Vc and X also depend implicitly on the ratio of matter density of isocurvature modes to that of 
the total matter at late time, i.e. whereas the various bispectrum contributions above 

do not. If 0.0/Vtm 1, it is possible to realise scenario where Vc Vn yet having 

7-)(cross2) 7-)(self2) i 

unsuppressed. 


3 Axion with gravitationally induced interactions 


Having discussed concisely the generic correlation structure between the bispectrum of the 
curvature perturbation B-ji and the cross-correlation power spectrum Vq in a generic two-held 
system, we now explicitly illustrate this generic behaviour using a concrete model. The model 
Lagrangian in the matter sector is m 


^ = \/^ 



- f^nf( 0 ) - Kxion(<h) - $) 


( 3 ) 


where 0 is the inhaton held and Hinf((/>) is its potential. We do not envisage any specihc form 
of Idnf(</>), as long as it can successfully support an inhationary epoch. The complex held $ is 
a PQ type held HU, which can be decomposed into real radial and angular components, r and 
9 respectively, as 


$ 


x/2’ 


( 4 ) 


The corresponding potential is the standard symmetry-breaking type, 14xion($) = A (1$^ — /f/2)^, 
so that the radial held is minimised at tq = fa with fa being the symmetry breaking scale, which 
is usually taken to be smaller than mpi. Assuming the radial held completely settles down at 
the minimum, we may identify the axion as the angular held a = faO. 

For concreteness, we consider a toy but generic dimension-5 interaction that can be gravi¬ 
tationally induced [12], 

= 9- -hh.c.. (5) 

mpi 

Additionally, overall we may add a cosmological constant to make the entire potential nearly 
vanishing at the global minimum. We shall associate the inhaton and axion huctuations with 
the curvature and isocurvature perturbations respectively in the standard manner as [13] 




FT , , 29 

——<p and X = — 

4>„ Do 


( 6 ) 


Here the subscript 0 denotes the homogeneous background held values. We have implicitly 
assumed Oe/Om ~ 1 as we are interested in the case where the isocurvature perturbations are 
large enough to be potentially observed in future experiments. 

Given the Lagrangian, we can straightly compute the correlation functions of the curvature 
and isocurvature perturbations using the in-in formalism JQ] [m [15] • The leading order contri¬ 
butions to the interaction Hamiltonian come from the gravitationally induced interaction term 
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in the matter sector, which at quadratic and cubic orders are 


H 2 = —a^a / d^x(j)6 with a = 2g —^sin(46*o), 


mpi 


^(seif 2 ) _ ^ 3 ^ I ^ -g^rUinUeo), 

J 3 mpi 

^cross 2 ) ^ f ^ 3 ^ = 4^ (.Os( 46 *o) • 

J "iPl 


_ 16 00 4 • 


H. 


(7) 

( 8 ) 
(9) 


The power spectra of the curvature and isocurvature perturbations, and their cross-correlation 
are 

n 2 


hV [h^ ^ 




7^x(A:) = - 


hV / 2 ^ ^ ^ 


27ry VVo 


2.,-3/2i^ (-fcr)'/'-"" 
r(3/2) ^ ^ 

2 ^ 9-312 r(^e) ('_^^^3/2-i.e 


r(3/2) 


Pc (A:) = 


7ra 


2rniy2 


A^/Vn{k)Vx{k ), 


( 10 ) 

( 11 ) 

( 12 ) 


where the indices of Hankel functions are dehned as = ^9/4 — with being the 

mass of 0 and similarly for ug. Here r is the conformal time dehned as dt = adr, whereas 
H and 0o are evaluated at horizon-crossing. Current observations suggest that the power 
spectra of isocurvature perturbations and the cross-correlation are much smaller than that of 
curvature perturbation on CMB scales, i.e. Vn 3> VcyVx, see the following. Here A is an 


integral involving the Hankel functions, dehned as ^ = 3? i dx/xH, 




x)h!)1\x) 


If we 


naively perform this integral to leading order in 3/2 and vg ^ 3/2 with a lower cutoh Xc, 
A ~ logXc so it is logarithmically divergent. But taking into account the asymptotic behaviour 
of the curvature perturbation at later times, we may evaluate A at an arbitrary value of Xc as 
long as exp [—1/C>(e)] < Xc < 1 [IS], giving A ~ —0.45. Details of the calculation are given in 
Appendix A. 

Next we consider the primordial bispectrum of the curvature perturbation that arises 
from the cubic interaction Hamiltonian. There are two contributions: one from the axion cubic 
self-interaction ~ 6^ and another from the cubic cross-interaction ^ Q‘^fp_ yyg 

therefore expect the resulting bispectra scale as for the former interaction and for 

the latter one. After some calculations similar as in ra. we hnd the resulting bispectra of the 
curvature perturbation are of the form, to leading order in the limit k, 3/2 and vg ^3/2, 




kl + kl + kl elP 

{hk2k^f ^ 3^3 ^4 ’ 

kl + kl + kl 2 /^ 167r3 ggy 
{kik2k^f ^3A^H^' 


(13) 

(14) 


Details of the calculation can be found in Appendix C. The above results of the two bispectrum 
contributions, Vc and Vxi are evaluated on relevant scales when the gravitationally induced 
interaction ([5|) becomes important, which do not necessarily correspond to CMB scales, but 
some smaller scales. 
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The amplitude of the bispectrum is usually measured in terms of the non-linear parameter 
/nl [IH]- The current constraints on /nl by Planck are = 0.8 ± 5.0, = —4 ± 43 

and = —26 ± 21 at 68 % conhdence level [2l] for k = 0.05Mpc. Having in mind /nl, we 

define a dimensionless shape function as 


r n 1 1 \ kik2k3 {kik2kzYBTz 


3 k\ + k2 + k^ {2 'kYV- 


n 


The corresponding shape functions for the bispectra and are then 

r(seif2)/T 1 k ) — and fh™ss2)/, u u \ _ 10 'Pc 

/nl and /^l (/ci, fcs,/cg) - “ 32 .^ 3/2 ’ 

' n 


97r^3 

or purely in terms of the model parameters as 


/, 


(self 2 ) 


NL 


H 


(15) 


(16) 


(roffO 


Again (fT 6 l) might suggest |/nl**^^^I < since Pc < P-r- This is not necessarily the 

case though as both of the non-linear parameters depend upon the coefficients 13 and 7 , which 
implicitly depend on the axion model parameters such as the misalignment angle Oq, as seen 
in flT7)) . It is possible to tune the axion model parameters such that |/nl^^^^| > |/nl°^^^^I while 
having Vc small. 

In the following we give an estimate about the size of contributions to the non-linear pa¬ 
rameter /nl on CMB scales in simple inflationary models. Assuming simple chaotic inflation 
with a quadratic potential, we have (po = 15mpi and the Hubble parameter at horizon-exit 
H PS From (ITTI) one then deduce 

/7ross2) _lg^ 4 gin 2 ( 4 ^o)cos( 40 o)(^) , (18) 


and I/nl^^^V/nl°^^^^I ~ fi'tan(46*o) sin(46*o)(</o/’^Pi)(^"o/-f^)^- Note that tq > H here since we are 
considering scenario where spontaneous symmetry breaking happens at higher energy scales 
than H. Therefore we expect I/nl^^^^I > I/nl°'^'^^^I simple example unless the coupling 

is very small, i.e. g < 10“^. In theory the value of /nl^^^^ could span a wide range depending 
on the combination of tq and g. For instance, assuming the symmetry breaking scale is around 
the GUT scale such that tq ~ 10“^mpi, then we find /nl^^^^ ~ 1 for ~ C>(10“^). 

From ffT 6 |) . we can also see that the resulting non-linear parameters are shape independent 
as opposed to non-canonical single field models such as K-inflation, where the non-Gaussianity 
induced due to non-linear interactions at horizon-crossing (or equivalently reduced sound speed 
Cs) typically peaks in the equilateral shape [19]. This is because the forms of non-linear inter¬ 
actions are completely different between our case and the single field case. An another example 
where non-Gaussianity induced by non-linear interactions at horizon-crossing does not peak 
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exactly in the equilateral shape is quasi-single held inhation [20]. In general models, however, 
the resulting non-linear parameters are not necessarily shape independent as we discussed in 
the last section. For our model, since /nl is shape independent, the tightest observational 
constraint comes from above if the gravitationally induced interaction becomes relevant 

on CMB scales. 

It is also convenient to egress Vc in terms of fractions of the isocurvature perturbation and 
cross-correlation, dehned aqj 


h = 


'Pn 


and Pc = 


Vc 


VVtiVx 


(19) 


For generic cold dark matter isocurvature perturbation, at k = 0.002Mpc“^ the upper bound 
from Planck TT, TE, EE -|- lowP -|- WP is Px ^ 0.021 [3|. On the other hand, the constraint on 
Pc from Planck TT, TE, EE -|- lowP -|- WP is —0.07 ^ Pc ^ 0.21 [3]. Recent study has shown 
\Pc\ needs to be of 0(0.1) for forthcoming CMB experiments to be sensitive to the isocurvature 
cross-correlation [10]. With the dehnitions flT^ . we can see how the observables are correlated 
in general 

/iJl ~and /ir"’~■ (20) 


3.1 Consistency checks 


There are implicit non-trivial constraints for the perturbation theory to remain valid, which 
also constrain the conditions for realising a large /nl given by (|T71) . In order to treat the 
transfer vertex a small perturbation, the correction to the curvature perturbation power 

spectrum Vn due to must be subleading compared to the leading piece flTOl) . Explicitly, the 
correction is given by 

APK(fc) = (-^) , (21) 

where V is an integral involving the Hankel functions dehned as 


V = ^ 


dxi 

Xi 




dx 




VQ 


'Xi 


X2 




( 22 ) 

with Xc being the lower cutoff, so that approximately "D ~ Thus for perturbative calcula¬ 
tions to be valid, from fl21l) we require 


7ia 


2rnFr2 


< 


0 ( 1 ). 


(23) 


Notice that this also means Vc < VV-rVx- For detailed derivations, see Appendix B. 

From (1T7|) . the non-linear parameter given in terms of the model parameter scales as 

oPp/ir^H^). Since the coefficients a and P are in fact related by /3 = — (40o/3)q; in our model, 
by using fl23l) and taking tq = fa-, we can see is as large as 0{1) only if 


pom-px 

p 

J a 

^Note this definition of is slightly different from /3iso used in [3]. They are related by /3i = diso/ (1 — diso). 
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This can be achieved by setting 0o ~ 0{mpi) and fa 0{mpi). Similarly, the coupling 7 
is related to the mass of the PQ held me- The massless axion limit, i.e. mg/H^ —)■ 0, thus 
corresponds to 


<^o7 

rlH^ 


< 1 . 


(25) 


From flTTI) . the non-linear parameter given in terms of the model parameter scales as 

/(rlH^). As a result, we can see ~ C>(1) only if 


mp\ 

00 


> 1 , 


(26) 


which is possible only if 0o is sub-Planckian. We can see the two conditions fl2T)) and (l26D are 
mutually exclusive in general. We therefore conclude only one of them can be made large and 
the resulting non-linear parameter /nl either scales as or Vq^/Vti- These results appear 
mainly due to the fact that the quadratic and cubic interaction Hamiltonians all arise from the 
same gravitationally induced interaction term (15|) . 

Our results also apply to more general cases where the gravitationally induced interaction 
between the inhaton 0 and the PQ held <F is of the form 

• ( 27 ) 

Uipi 

For m > 2, the cross-interaction 0^6* is also present, and contributes to /nl which is directly 
correlated with (3c 

~ (3c . (28) 

Whether or not this contribution can be comparable to and/or depends upon 0o, 

which is model-dependent. The relative size between and depends on the model 

of inhation with /nl°^^^~ ’^pi/ 0 O 5 whereas for the ratio go as ~ tan(?7,6*o). 


4 Conclusions 

In this article, we have studied the correlations between primordial non-Gaussianity and isocur¬ 
vature perturbation due to horizon-crossing interactions on general ground. In the presence of 
explicit couplings between the inhaton and another light held that later produces isocurvature 
perturbation, the bispectrum of the curvature perturbation receives contributions from the 
conversion of isocurvature perturbation. These contributions can be written in terms of the 
cross-correlation power spectrum. Taking the gravitationally induced coupling of an axion-like 
held as an explicit example, we have shown that the primordial bispectrum correlates with the 
cross-correlation power spectrum as B-ji Vc or Btz 'PcV'P'R, depending on the inhation 
model. If we do observe primordial non-Gaussianity on smaller than GMB scales in future 
experiments, we might be able to say something about or even put constraints on inhation 
models where the inhaton interacts with axion-like helds in the early universe. 

Besides the standard cold inhation scenario we have considered here, isocurvature perturba¬ 
tions can also arise naturally in warm inhation models via warm baryogenesis [22]. In such sce¬ 
nario, one should consider an inhaton-huid system instead. We hope to investigate the possible 
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correlations between the resulting isocurvature perturbation and primordial non-Gaussianity 
in such scenario in the future. 
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Appendix 


A. Computation of the two-point statistics 

Given the model Lagrangian ([3]), we can work out the corresponding Hamiltonian at appropriate 
order. At quadratic order, assuming slow-roll and the off-diagonal term of the mass matrix is 
small, we can take the free held Hamiltonian density, , and the interaction Hamiltonian 
density, 




H(,,e=-r + -wr + 


+ - 8a^g 


mpi 


cos(46'o)0^ 


7^2 = ~ 2 a^ 5 f—sin(46*0)06* = — 


mpi 


(29) 


where a is dehned as in ([71). The equations of motion of the free helds 0 and 6 follow from 

76 free- 

To compute the two-point correlation functions, we hrst promote the helds to operators and 
decompose them into annihilation and creation operators in Fourier space 


0k hktfk -|- , 

=bkVk + , (30) 


where Ok and fek satisfying the commutation relations 


4k, at = ik, it = (2ir)“(5<’>(k - k'), 


( 31 ) 


and zero otherwise. The mode functions Uk and satisfy the following equations 











where 




= 


and ml = —8g-^^rl cos(46*o). 


90^ 


mpi 


(33) 


Assuming the last terms in fl32|) are small such that both helds behave as massless helds up 
to slow-roll approximation, we hnd the solutions of the mode functions 


Mk(r) = — i exp 


2 


Wk(r) =-exp 

ro 


t {^9 + 


l\ TT 

7 2 

l\ 71 




2 / 2 




(34) 

(36) 


with — and similarly for ns, for 0 < n, < 3/2. Here is the Hankel 

function of the hrst kind. 

The correlation functions can then be computed using the in-in formalism |T3]. Using the 
commutator form, the expectation values of some operators 0{t) at the time t is given by 


o{t)) = (o 


Oit) 


0)+i dti ( 0 


f 




+ i dfi / df2 ( 0 




• £* 




0) + O{H^) 


(36) 


where |0) is the interaction vacuum at some initial time f*. The resulting power spectra ffTOj) 
flT^ can be found by inserting the interaction fl2^ into fl36l) for O = |0ki0k2) ^ki^k 2 ) 0ki^k2 
For instance, derivations of Vc and A can be found in |10j . 


B. Computation of the correction to the leading order Vjz 

In section 13.11 we argue that the condition fl23l) have to be satished for consistency in order to 
treat ^ small perturbation compared to 'Hfree- As stated in the section, this comes from 

the fact that the non-zero contribution to the inflaton power spectrum, or in other words, V-ji, 
coming from should be subleading compared to the free-held contribution flTOl) . 

To see, we compute the correction to Vtz due to using the commutator form of the 
in-in formalism (|36|) in a similar fashion as in UTl. Labeling this correction by APt^, this is of 
order ^[('H^)^] und is given by 


AVnik) = / dfi / df2(0 


-23fJ 


f*i 


dfi / df2 ( 0 


.Jt, 




(37) 


Assuming inhationary background such that a ~ {—1/Ht) and H th const, and taking the 
large wavelength limit after horizon-exit, we hnally arrive fl2Tll 


AV^^{k) 


|uk( 0 )|^ 




ira 


27r2 V 2rniL2 


V 


ira 

2rnH^ 


2 

VVn{k) , 


(38) 
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where T) is an integral involving the Hankel functions dehned as 

dxi . r dX 2 


V = ^ 


lo Xi 




In the massless limit where z/^, uq 3/2, V reduces to 


V = 


71^ 


3fJ 


dxi 


I 

'0 


[(1 — ixi)e^^^ — C.C.] (1 — fa:i)e 


fxi 


IXi 


X2 


<’(^2X’(i2) \ , (39) 


fxi 


dx2 

^2 


(1 + fa; 2 ) e 


2 —2ix2 


(40) 


and is logarithmically divergent in the IR limit. Here c.c. stands for complex conjugate. To 
regularsie this, we introduce a lower cut-off Xc in the above integral. Then we can see P 
approximately scale as 


V 


(log Xcf 
2 ! 




(41) 


For the correction to be small such that AV'iz{k)/VTi{k) <C 0{1), we therefore generically 
require | 7 rQ!/ 2 roi?^| 0{1), i.e. the condition 


C. Computation of the bispectrum 

The resulting primordial bispectrum can be computed in a similar fashion as the two-point 
functions. Assuming slow-roll, at cubic order, the corresponding interaction Hamiltonian 
dominated by the contribution from the gravitationally induced interaction ([5l) and is given by 

,3 ^ / Sxn, - , .3 4 . 3 ) / H H H 


X 


16 . . , 4 . , . 

-^00 cos (46*0 ) 6 'ki 6 ^ 1,2 6 'k 3 - - sin (46*0) 0ki 6*142 ^'kg + perm. 


(42) 

where perm, denotes permutation over momenta kj of the second term in the square bracket. 
Splitting in terms of the form of the interaction, we therefore arrive at ([HI) and ([9]). 

Cl. Contribution from 

Expressed in the nested commutator form of the in-in formalism, the leading contribution to the 
fully connected inflaton three-point function (0^) from the axion self-interaction Hamiltonian 
^^seif 2 ) cubic order in and is given by 




3\(self2) 


rt rti pt2 rti 

= dti dt2 / dts / dt4 

J to j to J to j to 


H\u), H\h), H^h), 


(43) 

where corresponds to an interaction Hamiltonian, with one being and the others 

being The overall contribution is therefore ini 


(^3)(..in) _j 2^3^ |«k,(0)'Ui,3(0)«k3(0)] 3} 


(*Tl 


f»T2 


rrz 


cLti 


dro 


dr?. 


'2 / U /3 / 

— 00 J —00 J —00 


dT^Wa^iji) 


2=1 


X (A + R + C)] (27r)3(5(^)(^k0 -|- 5 perm., 


(44) 
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where the terms A, B and C are 


A = - [Mki(ri) - c.c.][nki(ri)nki(T-2) - c.c.][nk3( 7 - 2 X 3 (r4)Mk3(r4) “ , 

(45) 

B = - [Mki(ri) - c.c.][Mk 2 (T- 2 ) - c.c.]X^(ri)n^^(r 2 )nki(r 3 )wk 2 (r 3 ) - c.c.]nk 3 (7-3X3 (7-4)43(^4) ’ 

(46) 

C =[Mki(ri) - C.C.][Mk 2 ( 7 - 2 ) - C.C.][Mk 3 ( 7 - 3 ) “ C.C.]4^ (ri)^^ (^ 2 X 3 ( 7 - 3 )Wki ( 7 - 4 )Wk 2 ( 7 - 4 )^^k 3 ( 7 - 4 ) , 

(47) 


corresponding to replacing H^{t 2 ), H\t^) or with the axion self-interaction Hamilto¬ 
nian respectively. Note that replacing H\ti) with gives zero contribution since 

oc which does not contract with any of the external legs 0/(t). Here we also have 
used the fact that Mk(0) is real. Plugging in the solutions for the mode functions u and n, 
fl45l) - (147|) written in terms of the Hankel functions are 


A = -W {-nf (-nf (-uf S 3 lHllH-km)HilH-hT2)] 


xS 


(48) 


B (-nf (-T 2 f (-nf » Hl'H-hn) S 


l'cl> 




(49) 


C = -iW(-ri)y-T2)y-r3)y-T4)®''"S Hi'Ji-hn) S Hil\-hT2) S 








x//®(-A^iT,)ff®(-A^3T2)ff®(-A^3T3)ff<;>(-;;ir4)i/«(-'=2r4)i/«(-*:3r4), 


(50) 


where we have dehned W = 8/rQ{H^/7^/2)^. Assuming TZ is sourced purely by the inflaton 
fluctuations, the resulting curvature perturbation bispectrum is simply 


(7^(kl)7^(k2)7^(k3)) 



(0(ki)0(k2)0(k3)) , 


(51) 


where the corresponding bispectrum k 2 , k^) and the dimensionless shape function 

^ 2 , ^ 3 ) as dehned in flT^ are given by 


(7^(kl)7^(k2)7^(k3)) =(27r)3<5(3) ( J^k, ) k2, k^) 


(2.)’5<3- (^k.) (4) (62) 


{kik2k‘iY 


Here H and 4 are evaluated at horizon exit. Factorising f^^^‘^\ki, /c 2 , k^) in terms of the shape 
dependent and independent parts, we have 




Xself2)s(self2)(A;i,fc2,fc3) 


[XYZfB 

X3 + y3 + Z3 ■ 


(53) 
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For convenience, we have introduced a normalisation wavenumber K and have dehned Xi = 
—Kxi, X = ki/K, X = k 2 /K and Z = fcs/JF in the above expression. In the massless limit 
^ 3/2, jg 

(self2) _ ( 9-3/2^9/2 


X' 


H 


U%v 


and the dimensionless shape function /c 2 , k^) is 


=3? 


r°°dXi f°°dX2 


Xi 


fxi 


X2 


•X2 


dX3 

X3 


dx. 


IXZ 


X 4 




X 


(x^f^ a [H>^J(Xx,)Hlf{Xx 2 )] S (Zx,) (ZX 4 ) Hj^ (Zx,) 

x//(;)(Kx2)ff® (yis) i/M (Via) 


+ {x,f ^-» (yijJ a {Yx^) hW(xx,)hi.I> (yi,)] 

xHil>(Zx3)Hif(Zx,)HS>(Zx,) 


i^(f) 


-i(xif^X XJ](Yx,,) » J/«(Za: 3 ) 




X <>(Jf3:iX> (ya:^) <> (Za:^) Hm(Xx,)Hm (Yx,) (Zx^)] } + 5 


perm. 


(54) 


(55) 


Here ^ 2 , ^( 3 ) diverges in the IR in the massless limit. Introducing a lower cut-off Xc as 

before, one hnds 





TT ^/"^{logXcf 


X3 + F3 + Z3 
{XYZf!'^ ’ 


(56) 


to leading order. In the following, we briefly summarise how to compute the above result. 

To compute the leading order term in the IR limit, we take the following asymptotic form 
of the Hankel function in the a; 1 limit 


( 1 ), , .23/2r(3/2) _ 3/2 .2-i/2r(3/2) 1/2 

n Jx) X - I - —!—^X - I - , ; ' ' X^'^ 

3/2^ ^ TT 7r(l/2) 

/ ^.cos( 37 r/ 2 )r(- 3 / 2 ) 1 \ 

V * 23/% ^ 23/2r(5/2) J 


x3/2 + ... 


(57) 


whereas that of H^j^x) is the complex conjugate of (|57)) . Here we can see the real and the 
imaginary parts scale as 0{x^^‘^) and 0{x~^^‘^) respectively at leading order. Using this, we can 
see how the integrand in fl55l) scales with Xc- For the second line in fl55l) . at leading order in 
the IR limit, it scales as 


i=l 


HilHXx,) 


—^ H(^3 ^ 

2=1 


1 

23/2F(5/2) 


(Xxi)3/2 + ... 


(58) 
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whereas for the third line, the leading order piece is 


- 3 (ZX 2 ) <> {Zxi) Hlf (Zxi) 


-)■ 


r(3/2) 

_r( 5 / 2 )ir 


TT 




XiJ 


( 59 ) 


Here we have taken {Xxi), hII^ {ZX 2 ) and either {ZX 4 ) or (Zx^) to be imaginary. 
For the fonrth line in (|5^ . it must be real in order to keep the whole integrand real. The 
leading order term scales as 


nil' (Yx2)Hlf (Yx3)Hil>(Yx3)^2 


r(3/2) 1 

■23/2F(3/2)' 

_F(5/2)7r_ 

TT 


{Yx-2)-Y'‘ + ... 


(60) 


Here we have taken iXX 2 ) and either (Yxj) or {YX 3 ) to be imaginary. Together 

the above leading order pieces in the three lines therefore give the following integrand 


^=1 


r(3/2) 

_r(5/2)7r_ 


■r(3/2)2^/2- 

TT 



(61) 


and the corresponding contribution to g^^les as (log Xc)'^ in the IR limit. Taking 

other choices of the leading order pieces in the Hankel functions give subleading terms. For 
instance, taking hI) 1 \Xxi) to be real and Hul\Xx 2 ) to be imaginary instead in the second line, 
together with the above leading order terms in the hrst and third lines, the resulting integral 
scales at most as (—logXc)^. 

Similarly, for the hfth line in fl55l) . the leading order piece is 


S fcp (Kx2)j a [H'Y>(Xx,)Hm {Yx2) H<1\Xx3)HW (Kxs)] 




r(3/2) 

_F(5/2)7r 


23/2f(3/2) 


TT 


(Xx,) 


- 3/2 


Y \ 

f 


Y \ 

x) 


(62) 


Here we have taken Hu^\Xxi), {YX 2 ) and either H^g[Xx^) or iYx^) to be imaginary. 
In order for the overall contribution together with the second and hfth lines to the integral to 
be real, the leading order piece in the sixth line must be real and is given by 


//<;> (Z13) //<;) (Z14) (Z14) ^ 


r(3/2) 1 

'2Y2r(3/2)' 

_F(5/2)7r_ 

TT 


(Zx3)-Y^ + ... 


(63) 


Here we have taken (Zxs) and either (ZX 4 ) or (ZX 4 ) to be imaginary. To- 
gether the second, hfth and sixth lines in (l55l) give the following leading order contribution 
to ^ 2 ,^ 3 ) in the IR limit 


i=l 


1 

r(3/2) 1 

4 

'2Y^r(3/2)' 

Y ^ y /2 X y y/ 2 - 


_F(5/2)7r_ 


TT 

\Yz} ^\Xz} 


(64) 
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Again similarly, the leading order term in the seventh line in fl55l) is 


z 3? {Yx^)] 3? [i/W {Zx^)] ^ -i{x^f/^ 


l^cj) 


_23/2r(5/2) 


{Yx2f/HZX3f^‘' + ... 


(65) 

whereas for the eighth line, the leading order term must be imaginary in order to keep the 
overall contribution real and thus is given by 


(Yx,) <> (Zx,) XX(Xx,)Hm (Yx^) »<« (Zx^) 

^ y/2 / y \ 3/2 / y \ 3/2' 


I 


'23/21(3/2)' 

4 

r(3/2) 1 

TT 


_r(5/2)7r_ 


{XY Z)-^/\xiX2X3)-^/^ 


+ 


( 66 ) 


Here we have taken either Hj)l\Xx 4 ), iYX/^) or (ZX 4 ) to be real. Together the second, 

seventh and eighth lines give the following leading order contribution to /c 2 , /cs) in the 

IR limit 


H' 

i=l 


Xi 




r(3/2) 


r(5/2)7r TT 


23/2r(3/2) 


X \ 3/2 / V \ 3/2 / 7 \ 3/2 


(67) 


After factorising out the common factor 11^=1 (^*) ^ ([6T]), fIMl) and (1671) . we can see the 

integral goes as 


4 

dxi / dx2 / dxa / dx^ 
Jxi J X2 j X3 


xd ^ ~ 


(log XcY 
4! 


( 68 ) 


assuming the integral vanishes in the UV limit —)■ 00 . Thus dually, to leading order in the 

IR limit, /c 2 , /cs) is given by 


^{XYZ)-^/^ 


IR 


3/2 

r(3/2) 1 

4 

'23/21(3/2)' 


_r(5/2)7r_ 


TT 


(logXc 

4! 


X [dX^ + 2{X^ + y3) + (^3 + y3 y z^)] + 5 perm.. 


(69) 


Summing over all permutations, it is not difficult to see /c 2 , ks) scales as (X3 + y3 _|_ 

Z^)/{XY and thus f^hiki, k 2 , is approximately shape independent. Using the 

results of the cross-correlation spectrum Vc and expressing fT<fh{ki, /c 2 , in terms of Vc, 

we therefore arrive at (HE!). 

C2. Contribution from 

The corresponding leading order contribution to the inflaton bispectrum {(pY from the axion 
cross-interaction ffh™ss 2 ) quadratic order in Expressed in commutator form, this 

contribution is given by 


’>'"“■")=/ *1 / dh I d(3([i/'((3), [i/'((2),[i7'((l),.).,(()=]]]> , 


rh 


f*t2 


(70) 


'io 


'to 


'to 
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where one of the is given by 7 ;fh™ss 2 ) Expanding this, we get 


(03)^"°"s2) ^ _ 4 q, 2 ^ j j (27r)^5® (ki + ka + kg) 


.2=1 


x3? 


—t 


dri / dra / drg a^(ri) (D + E + F) 


' —OO J—OQ J—OO 


2=1 


+ 5 perm., (71) 


where D, E and F correspond to contributions coming from replacing 77/(fa) and 

respectively with f/j^^cross. 


D = [ukAn) - c.c.][vlAri)vkA'r2)uk2{r2) - c.c.]vkAE)vlA'rAlulA'rAl , (72) 

E = [Mfci(ri) - c.c.][nfc^(ri)nfci(ra)nfc2(ra) - c.c.]nfc3(ra)nfc3(r3)Mfc3(r3), (73) 

F = [Mfci(ri) - c.c.][Mfc2(ra) - c.c]vkA'ri)vk2{r2)vlArAK2i.^^)KA'^^) ’ ('^4) 

Here we have again used the fact that Mfc(O) is real. Written in terms of the Hankel functions, 
they are given by 


E = - iN(-nf(-T2fH-T3f» [i/lp-fcin)] S [i/if 

-yH<^'’(-k3T2)HilH-k,T,)HilH-hT,), 


(75) 

(76) 

(77) 


where M = 4/ro(7/y^/2)^. The resulting the curvature bispectrum k 2 , ks) is sim¬ 
ply given by multiplying (171 II with —(/7/'^o)^- Similar to the case, we factorise the 

correspon(iing shape function k 2 , k^) into shape independent and dependent parts 

as follows 


fNLih, k2, A)3)(“°^^') = ^(--2)^(cross2)^^^^ 
In the massless limit ug, ^ 3/2, j^(cross 2 ) jg given by 


(XYZfl‘^ 

X3 + y3 + ^3 • 


(78) 


^(cross2) 



CK^7 ^ 

Hh^J V2J ’ 


(79) 
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whereas k 2 , k^) is 




=3? 


lo Xi 


'Xl 


dX2 

X2 


•X2 


dx3 

X3 


3f? 


HllKXx,) 


3/2 


3f? 


H<;f(Yxi)HW (Kx2) 7/W (Kx2) 


xi/m (Zx,) (Zx,) (Z 13 ) - ( 13 )=’-'^ S (Kxj) 

xi/P) (Zx2) i/W (Zxs) (Z13) - (13)=’-'^ S [i/W (Ki^)] H<^\Xxi)H!.f (Yx,) 

X XX{Xx,)Hl'J (Yx,) H'X (Zx,)]} + 5 perm. 


(80) 


Here we have again dehned Xi = —KTi, X = ki/K, Y = ^ 2 /-^ and Z = k^/K. Similar to 
the case, /c 2 , ^ 3 ) diverges in the IR in the massless limit. Introducing a lower 

cut-off Xc and following similar approach as in the case, one hnds 


S^^^^^^^\ki,k2,k3)^-{XYZ) 

IR 


3/2 

r(3/2) 1 

3 

■23/2r(3/2)' 

■(log Xcf 


_P(5/2)7r_ 


71 

3! 


- 2{X^ + Y^) - {X^ + + Z^)] + 5 perm.. 


(81) 


Summing over all permutations, it is not difficult to see ^ 2 , k^) also scales as -|- 

Y^ + Z^)/{XYZfl^ and thus /NL(fci, ^( 2 , is approximately shape independent. Rewrit¬ 

ing f-tih{ki, ^ 2 , /c 3 )h™s® 2 ) terms of Vc, we therefore arrive at flThl) . 
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